Quantum decoherence of a charge qubit in a spin-fermion model 
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We consider quantum decoherence in solid-state systems by studying the transverse dynamics of 
a single qubit interacting with a fermionic bath and driven by external pulses. Our interest is in 
investigating the extent to which the lost coherence can be restored by the application of external 
pulses to the qubit. We show that the qubit evolution under various pulse sequences can be mapped 
onto Keldysh path integrals. This approach allows a simple diagrammatic treatment of different 
bath excitation processes contributing to qubit decoherence. We apply this theory to the evolution 
of the qubit coupled to the Andreev fluctuator bath in the context of widely studied superconducting 
qubits. We show that charge fluctuations within the Andreev-fluctuator model lead to a 1/f noise 
spectrum with a characteristic temperature depedence. We discuss the strategy for suppression of 
decoherence by the application of higher-order (beyond spin echo) pulse sequences. 
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I. INTRODUCTION 

The loss of coherence of a quantum two-level system 
(quantum bit) is caused by its unavoidable coupling to 
the surrounding environment. For solid-state qubits, the 
decoherence process can be quite fast due to coupling 
to a large number of internal degrees of freedom. Our 
understanding of quantum decoherence and methods for 
its suppression in a realistic solid-state environment is 
mainly confined to the cases of a qubit interacting with 
bosonioii 2 -^ and nuclear spin baths^^ the so-called (and 
extensively studied) spin-boson and spin-bath models, re- 
spectively. A less well understood, but very relevant case 
for solid-state quantum architectures is that of a qubit 
coupled to a fermionic bath^&^ i 10 ' 11 ! 12 ! 13 which dramat- 
ically differs from the previous examples. In this paper 
we study quantum decoherence in the context of a su- 
perconducting charge qubi t 14 ' 15 i 16 i 17 i 18 ] 19 i 20 interacting 
with the non-trivial bath of Andreev fluctuatorsj 21 ' 22 i 23 
This problem is a paradigmatic spin-fermion decoher- 
ence problem and applies to many situations involving 
the quantum coupling of a qubit ("spin") to a general 
fermionic environment. Using a many-body Keldysh 
path integral approach ) 9 ' 10 i 24 i 25 i 26 we obtain a quantum- 
mechanical description of the qubit evolution under pulse 
sequences aimed at prolonging the coherence of the sys- 
tem. The simplest case of decoherence under pulses 
is the spin echo dephasing experiment, which has been 
shown to extend the coherence time of solid state (su- 
perconducting) qubits ; 15 ' 17 ! 27 by essentially eliminating 
the quasi-static shifts of qubit energy splitting (inho- 
mogeneous broadening) due to the slow environmental 
fluctuations. However, sequences involving more pulses, 
for example, CPMG 2 ^ and Uhrig' o 2 ' 29 sequences, are ex- 
pected to lead to a further increase of the coherence 
time . 2 - 29 - 30 - 31 - 32 - 33 - 34 - 35 ' 36 - 37 

In this paper, we consider an experimentally relevant 
example - a superconducting qubit coupled to fluctuat- 
ing background charges ) 18 ! 19 e.g. electrons residing on 



Anderson-impurity sites. Due to a large on-site Coulomb 
repulsion forbidding double occupancy, this example rep- 
resents a non-trivial interacting bath. The dynamics of 
the charge fluctuations on the impurity sites is deter- 
mined by the hybridization of impurity levels with the 
quasiparticle band of the superconductor. To the low- 
est order in tunneling at the superconductor/insulator 
interface, the hybridization of the impurity levels can 
be described by a correlated tunneling events of two 
electrons with opposite spin to/from the superconduc- 
tor. We show that in the small background-charge den- 
sity limit, these fluctuations lead to a 1/f spectral den- 
sity of noise. Using these results, we finally obtain the 
quantum-mechanical description of the qubit evolution 
driven by external pulses, and discuss optimal strategy 
for the suppression of the decoherence with designed com- 
posite pulse sequences. 

The paper is organized as follows: In Sec.|IT]we provide 
a general derivation of the qubit evolution with pulses 
and map the calculation of decoherence function onto 
Keldysh path integral formalism. In Sees. IIIII and IIVI 
we introduce Andreev fluctuator bath and derive spec- 
tral density of noise for this model. Finally, in Sec. |V| 
we discuss the influence of pulse sequences on the qubit 
decoherence. 



II. GENERAL THEORY FOR QUBIT 
EVOLUTION 



The transverse dynamics of a qubit interacting with its 
environment is determined by the following Hamiltonian 



H 



E 



a z V + H B ■ 



(1) 



Here the environment is represented by a fermionic bath 
Hb , and the qubit is coupled to the environment through 
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the density fluctuation operator: 



V = 



(2) 



la 



This model corresponds precisely to the coupling of a 
superconducting charge qubit to the density fluctuations 
on the impurities in the substrate. Here c; CT and c[ a are 
the fermionic annihilation and creation operators at Z-th. 
site with spin a, and vi and (ni a ) are, respectively, the 
strength of the coupling and average occupation of Z-th 
impurity, i.e. (ni a ) — (cj^Qo-). Equations (p} and ([2|) 
define our spin-fermion model. 

We study the evolution of the qubit in contact with 
a fermionic bath assuming the qubit energy relaxation 
time T\ to be much longer than the quantum dephasing 
time Ti (thus only a z coupling is present in the Hamil- 
tonian). Qubit decoherence under the influence of the 
environment is given by the off-diagonal matrix elements 
of the qubit 's reduced density matrix, and for the free 
evolution of the qubit we get (h = 1) 

P+ -(t) = <+|Tr B {p(*)}|-) = p + .(0)e- iBt W(t). (3) 

In the above p(t) is the density matrix of the whole sys- 
tem (qubit+bath), which is assumed factorizable at t=0, 
Ttb{...} is the trace with respect to the bath degrees of 
freedom, and W(t) is the decoherence function defined 

as 



W(t) 



= /e l(HE 



-V)t -i(H B +V)t 



(4) 



with the brackets representing the thermal average with 
respect to the bath Hamiltonian Hb, i.e. (...) = 
Tyb{pb---}- The time t always refers here to the total 
time of the evolution. 

In addition to the free evolution of the qubit 
(free induction decay), one is often interested in 
the dynamics of the system subject to external 7r- 
pulses^^^^i^i^i^i^i^i^l which could, in principle, 
prolong or restore quantum coherence. The 7r-pulses con- 
sidered here correspond to rotations of the qubit 's Bloch 
vector by angle n about, e.g., the x axis, and are short 
enough for the bath dynamics during the pulse duration 
to be negligible. Then, the evolution operator for qubit 
and bath is given by 



iHr 



iHr 



,-iHri 



(5) 



with n and r< being the number of applied pulses and 
time delays between the pulses, respectively, and the to- 
tal evolution time t = J^'i T *- One can see that the 
well-known Hahn spin echo (SE) sequence, for example, 
corresponds to a single pulse with t\ =ti =t/2. 

Using the fact that in the "pure dephasing" case under 
consideration, the qubit states |±) are the eigenstates 
of the Hamiltonian (JTJ, we can write the decoherence 
function under the influence of pulses as 



W n (t) = ( (C/i n) (t)l U™(t) 



t „ 



(«), 



(6) 



with the evolution operators (t) given by 
J7| n) (t)= e - l (^ B + 1> ) T "+ie _4(i * B_ ' t>)r "...e"' ( ^ B+pt>)Tl 



J7^" ) (t)= e - l (^ B - t> ) r "+ 1 e - l (^ B+1> ) r "...e- l( ^ B ~ !3 ^ )7 " 1 



(7) 



where p = (—1)" is the parity of the sequence. Then, 
the off-diagonal elements of the qubit density matrix are 
given by 



P+ 



_(t) - p Pl - p (0)e- i P E ^- T *+~ + P T "+^W n (t). (8) 



Here the phase factor is zero for all balanced sequences 
(for which the total times of evolution due to H+V and 
H — V are the same in Eq. Q). The evolution of the 
qubit under SE sequence, for example, acquires a simple 
form 



P s + E - (*) = p-+ (0) ( e m + * e m - i e - iH +h- iH - 



(9) 



with H± = H B ±V. 

Decoherence under pulses has been analyzed with 
methods specific to the spin-boson model^i and the spin 
bath model,— & or using operator algebra . 33 ! 34 ' 36 The lat- 
ter approach, although very general, does not allow for 
transparent understanding of physics of the bath. How- 
ever, the evaluation of W n (t) defined in Eq. (JBJ) can be 
mapped onto the evolution on the Keldysh contour, 38 
putting the calculation of decoherence into the frame- 
work of many-body theory. Similar formalism has been 
used to study full counting statistics of a general quan- 
tum mechanical variable and has proved to be quite 
convenientj 24 ' 25 The evolution operators fj± can be 
written as 



U { ±\t) =Texp 



-i / (H B ± f n {t')V)dt' 



(10) 



where T is the time ordering operator. The function 
f n (t') encodes a particular sequence, and is defined as 

n 

U(t') = P ^(-i) fc e(z fe+1 - t')e(f - t k ), (ii) 

where Q(t') is the Heaviside step function, tk with k — 
l..n are the times at which the pulses are applied, t — 0, 
and in+i = t. Thus, the product of operators inside the 
average in Eq. ((6|) corresponds to (reading from left to 
right) the time-ordered evolution from to t (with +V 
coupling), followed by the time anti-ordered evolution 
from t to (with — V coupling). We can then introduce 
the Keldysh contour C (see Fig. [5^,) together with the 
notion of contour-ordering of operators ! 26 ' 38 The qubit- 
bath coupling takes then two opposite signs on the up- 
per/lower branch of the contour: Vc=±V. While f n (t') 
is non-zero only for t' £ [0, t\, we can extend the limits 
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of time integration on both branches to [—00,00]. The 
evolution from t' = — 00 allows one to include the adia- 
batically turned-on interactions in H B (see. for example, 
Ref. [2y|), paving the way to the treatment of decoher- 



ence in an interacting fermionic bath. The final result 

is most compactly written as a functional integral with 

the Grassmann fields 4>i and -0; defined on the Keldysh 
contour&io^^s 



W n (t)=(T c exp 



-11 dt' 
'c 



H B +Vcfn(t') 



Z B 



T^iThl)iexpiiS B [^,fl>]-i / dt'^i (*')/» {t')[^a{t')M') ~{n la }]\ 



la 



(12) 




FIG. 1: (color online). Correlated tunneling of two electrons 
with opposite spins from the impurity sites in the insulator 
into the superconductor. An electron from the i-th impurity 
with energy below the gap A tunnels into superconductor, 
propagates over distances of the order of coherence length 
£ and recombines with another electron with opposite spin 
from j-th site into a Cooper-pair. The amplitude for such 
Andreev process decays exponentially with distance between 
the impurity sites Ay <x exp(— |rj — Tj\/it£), see Eq. (TH)) . 



where the integration is performed on the contour C 
shown in Fig. [2^, vi(t') — ±vi on the upper/lower branch 
of the contour, and the normalization constant is defined 
as the functional integral with V = 0. The bath action 
Sb = So + Si n t , and the functional integration with non- 
interacting Sq corresponds to averaging over an equilib- 
rium noninteracting density matrix at t' — —00. This 
formulation of the decoherence problem enables one to 
use techniques and approximations developed in many- 
body theory. It also allows for a transparent treatment of 
the physics of the bath while simply encoding the driving 
of the qubit in a single function of time f n (t'). 



III. ANDREEV FLUCTUATOR BATH 

In order to evaluate the functional integral (|12[). 
one needs to specify the bath Hamiltonian. Here, 
as an example, we consider a non-trivial bath of An- 
dreev fluctuatorS ) 21 i 22 ' 23 which describes the fluctua- 
tions of the occupation of impurities close to insula- 
tor/superconductor interface due to Andreev processes. 
This model takes into account coherent processes of cre- 
ation (destruction) of the Cooper pair in the supercon- 



ductor by correlated tunneling of two electrons from 
(to) different impurity sites in the insulator , 22 ' 39 see also 
Fig. [TJ In the limit when the superconducting gap en- 
ergy A is the largest relevant energy scale in the prob- 
lem (T, E,Ej,<g. A), the effective Hamiltonian for the 
Andreev fluctuator bath, after integrating out supercon- 
ducting degrees of freedom, is given by 



la I 



(13) 



Here, Si and U are the energy of a localized electron 
on l-th impurity (measured with respect to the Fermi 
energy £j? of the conduction electrons) and repulsive on- 
site interaction (assumed to be large enough to prevent 
double occupation of the sites), respectively. The matrix 
elements Ay, in the limit of low transparency barrier 
between the insulator and superconductor, are given by 



A 



ij 



A 



sin(p F |r ; - r-j 
p F \ri - rj\ 



|ri-rj|/7r£ 



(14) 



Here p F is the Fermi momentum, £ is the coherence 
length in a clean superconductor. The amplitude A = 
2ir 2 d 2 aN (Q)T Q 2 is determined by the tunneling matrix el- 
ement between the insulator and superconductor To, the 
normal density of states in the metal N(Q) — mpp/n 2 , 
the localization length under the barrier d and the size 
of the impurity wavefunction a.- 22 

Given the Hamiltonian (fT3")) , the action for the bath on 
the Keldysh contour can be written as 

13 I 

+A* lj i> n {t')$ n {t') + Ay^;(/)<M(*') ■ (15) 

Here we used the mean-field approximation for the An- 
derson impurity model assuming that the Kondo tem- 
perature Tk is smaller than the superconducting gap A, 
which is reasonable in the situation at hand when the 
impurities are located in the substrate and the tunnel- 
ing matrix element Tq coupling them to the states in the 
superconductor is small. The occupation probabilities 
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{ni a ) are obtained self-consistently using 

M = / £ n H")[G&>) ~ (16) 

see Ref. for more details. Performing a Keldysh 
rotation ; 26 ' 38 one can calculate the full Green's function 
G Ua {t,t') for the bath (see Fig. [3]) 



the amplitude Aij oscillates on the length scale of Pp 1 , 
the contribution of these off-diagonal terms to the self- 
energy is small. 

Using the above results, the action for the bath can be 
written in terms of the full Green's function Gu a (t,t'). 
Then, the decoherence function becomes 



G^(t,0=^(*.0-SH t r(*,0- 



(17) a) 



Here G u ^(t, t') is the bare Green's function, see Eq. (fl~5)) . 

and the self energy Y,u a (t, t') is calculated to second order 
in giving the components of the self-energy matrix 

= ^lA^Gflti^t), (18) 

= Ei^i 2 ^-^*'-*)- ( 19 ) 

In Eqs. (HHJ) and (TT5)) we have neglected the off-diagonal 
terms in the impurity indices, i.e. £y CT ~ bij^iia- Since 



i' = 



6) 



1 -- 



/i(*'<0) = 



V C = V 



V c = -V 



f=t 



/i(*'>0)=0 



f 



FIG. 2: a) Dependence of Vc(t) on time along the Keldysh 
contour, b) The plot of the function f n (t') for the Spin Echo 
sequence (n = 1) . 



W n (t)=exp[- X n(t)} = 



x exp 



Z B 

oo oo 



(20) 



E Mi / *=» E 



-oo — oo 



,6=1 



Here, Pi CT (t) corresponds to the fermion density opera- 



tor pia(t) 



; the fields 



rpW(t) and ip^ 2 \t) are given by the appropriate super- 
position of the fermionic fields on the upper and lower 
parts of the Keldysh contour, see Ref. [26|]. After per- 
forming the functional integral over the fermionic fields 
and expanding to second order in vi, one finds 

, ^JoJo 



+G^(t 1 ,t a )G^(t a ,t 1 )+G^(t 1 ,t a )G^(t 2 ,t 1 )J.(21) 

Equation (|21j) holds whenever the short-time expansion 
is valid. The long-time asymptote can be obtained by 
resumming the whole seriesi*^ 1 ^ 

By introducing the Fourier transform of the Green's 
functions, Eq. ([21]) can be formally recast as 



Xn(t) = 



duj F n (u)t) 
2^ w 2 



(22) 



Here F n (uit) — uj 2 \f n (u>)\ 2 /2 is a sequence-specific filter 
function, the role of which we discuss in Sec.[V] Thus, to 
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FIG. 3: (color online). Dyson's equation for the retarded 
Green's function G R in the Born a ppr oximation. Here we 
have adopted the convention of Ref. [38(]. The advanced and 
Keldysh Green's functions are obtained analogously resulting 
in Eqs. Ip))-([I51). 



second order in Vi we obtain Xn(t) having the same struc- 
ture as in the case of a qubit coupled to the spin-boson 
bath or classical noise,— ' 29 ' 37 ' 41 i.e. Xn(t) is the integral 
of the product of the environment-specific spectral den- 
sity of noise Sq(lo) and sequence-specific filter function 
F n (wt). The spectral density of quantum noise Sq(uj) in 
the spin-fermion problem is given by 



5 



S Q (U): 



la 



2^ 



(o-|)+g« (n+|) el (n-|) + o£ (o+|) g- (n-|)] . (23) 



In the frequency domain, the full Green's functions are 



1 



U-El- U{ni-a) - £^ H (u>) 

G&H = tanh (^) [g£» - G#»] 

where the self energy is defined as 

vA/R, a _ I A? 1 2 



oj + Ej + U (rtjcr) T ^ 



(24) 



(25) 



Equation (|23|) . defining the noise spectral density in the 
quantum-mechanical many-body language enables a di- 
rect calculation of decoherence in various situations, as 
we consider next. 



IV. SPECTRAL DENSITY OF NOISE DUE TO 
ANDREEV FLUCTUATORS 

In general, the solution for Sq(lu) with many Andreev 
fluctuators, can be obtained numerically by randomly 
generating the energies £/, and positions r; of the im- 
purities at the insulator/superconductor interface. The 
numerically obtained spectral density of noise Sq(cj) for 
50 fluctuators is shown in Fig. [U At low frequencies the 
noise power spectrum has 1// dependence. 

For u) and Aq much smaller than the typical impurity 
level spacing 6s and temperature T, the analytical solu- 
tion for the spectral density of noise (f2"3"]) is given by 



Sq(w)i 



1-tanh 



2T 



liable) / 9fi x 

w 2 + 4 7iff (e ;CT ) 2 ' [M) 



where Ei a = £i + U (n;,_ CT ), and ji a (itv) = ImE^sv) is 
the broadening of the impurity energy levels due to An- 
dreev processes. This broadening corresponds to the fluc- 
tuations of the impurity occupations changing the elec- 
trostatic environment of the qubit, and thus causing de- 
phasing. From Eq. ([26]). one can see that Sq(u) is given 
by a sum of Lorentzians with different widths, which un- 
der proper distribution of 7 ; gives rise to a 1// noise spec- 
trum (see below) . Given that the charge density fluctua- 
tions via Andreev processes involve two impurities with 
energies of the order of T, the probability to find two 
such impurities is proportional to (T/D) 2 with D being 
the impurity energy bandwidth, and thus, Sq(u>) oc T 2 
at low frequencies as seen exp er iment ally . 1 8 ' 1 9 

For 1/ / spectrum to arise from Eq. (j2"rJ]l . the distribu- 
tion of 7; has to be log-normal. In order to have such dis- 
tribution, the density of the charge traps has to be small, 



Sq(0) 

0.500 



0.100 
0.050 



0.010 
0.005 




0.1 



0.2 



0.5 



1.0 



2.0 



FIG. 4: (color online). Log-log plot of the noise spectral den- 
sity Sq (u>) for Andreev fluctuator model. The plot is obtained 
by randomly generating the positions r; £ [—5, 5]£ and ener- 
gies Ei £ [—1, 1]K of 50 impurity sites and then, numerically 
integrating Eq. ((23} • Here we assumed that Vi = v, on-site re- 
pulsion U — > 00, and the sites are occupied with equal number 
of electrons with spins up and down. We used pj, 1 — 10~ 2 £, 
A) = 0.1K and T = 0.1K. The dashed (blue) and dot-dashed 
(black) lines, shown for comparison, correspond to 1// and 
1// noise spectra. 



so that the dominant contribution to the self energy in 
Eq. (|25p comes from few pairs of impurity sites, which 
are selected from the sum because of the energy conser- 
vation and distance constraint. Then, the switching rate 
7/ oc exp(— 2|r; — rj |/7r£) for a certain j (see Eq. (TT^l). 
Since the distances between the charge traps are uni- 
formly distributed, the probability of finding a switching 
rate 7 is P( 7 ) oc l/ 7 , leading to 1// noise. In the op- 
posite limit of large density of charge traps, many sites 
j contribute to the sum in Eq. (|25p . and the switching 
rates 7/ self average and become approximately the same 
for all sites. Note that unlike in the theory of 1// charge 
noise produced by fluctuating two level systems (TLS) in 
the substrate with log-uniform distribution in the tunnel 
splitting^ the emergence of the 1// noise within An- 
dreev fluctuator model has a qualitatively different geo- 
metrical origin due to the exponential dependence of the 
rate 7; on the distance between different impurity sites. 
This finding of the geometric origin of the 1// noise in 
the Andreev fluctuator model is an important result of 
our work. 

We note that the model of charge traps with no on- 
site repulsion U = [2l[ does not lead to 1// noise be- 
cause in this case the self-energy is dominated by the two- 
electron tunneling from the same site. The contributions 
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to the self energy from Andreev processes involving other 
sites are exponentially smaller than the dominant term, 
and the distribution of the rates in Eq. |26|) is not log- 
normal. Therefore, we emphasize that the realistic model 
for 1// noise due to Andreev processes should include 
both spinful fcrmions (to correctly describe the dynam- 
ics of charge fluctuations), and large on-site repulsion (to 
prevent double-electron occupation). 

At high frequencies u> ^> Se, T, the spectral density 
Sq(uj) has resonances corresponding to the virtual pro- 
cesses of correlated two-electron tunneling from (to) the 
impurity sites in the insulator. These resonances, de- 
scribing manifestly quantum-mechanical processes, can 
be seen in Fig. [4] at high frequencies. Their contribu- 
tion to the decoherence of the qubit is suppressed by a 
factor F n (ujt)/uj 2 , see Eq. (j2"2")l . However, going beyond 
the pure dephasing model, Ti3> T2, considered here, one 
can show that correlated two-electron tunneling processes 
contribute to the energy relaxation of the qubit «2i 

V. THE INFLUENCE OF PULSES ON 
DECOHERENCE 

The time dependence of the decoherence function 
W n (t) under a pulse sequence is given by Eqs. (|2"0)) - 
(22), showing that the noise contribution is modulated 
by a filter function F n (u>t). For the free induction de- 
cay Fo(u>t) = 2 sin 2 [wt /2] , whereas for spin echo we 
have Fi(uit) = 8 sin 4 [wi/4] suppressing the low- frequency 
(w<4/t) part of Sq(lu). In general, higher-order pulse 
sequences act as more efficient high-pass filters of envi- 
ronmental noise, i.e. for n pulses applied in time t only 
frequencies Lu>2n/t contribute to x n (t)- Due to the for- 
mal analogy between Eq. (|22[) and the solution for the 
decoherence under classical Gaussian noise, the analysis 
given for the latter case in Ref. [37| also applies here as 
long as the time expansion is valid. The results rele- 
vant for the noise spectral density derived here can be 
summarized as follows. For Sq(lu) oc 1/w q , we obtain 
Xn(~t) K t 1+a /n a for all sequences applicable to the pure 
dephasing case, i.e. the CPMG sequence^ periodic dy- 
namical decoupling, 3 — concatenations of spin echO ) 33 i 34 
and Uhrig's sequenced Thus, sequences beyond spin echo 
should lead to a further increase in coherence time for 1 // 
spectral density of noisei 32 ' 35 ' 37 For noise spectrum with- 



out sharp ultra-violet cutoff, which is the case considered 
here, the CPMG sequence marginally outperforms other 
sequences! 29 ' 37 Furthermore, taking into account the sim- 
plicity of CPMG sequence (defined by t\ — r n+ i = t/2n 
and all the other Tj = t/n), we believe that it is a pre- 
ferred approach of noise suppression for the problem at 
hand. 3 ' We therefore propose that detailed experimental 
investigation of superconducting charge qubit dephasing 
behavior be carried out in order to test our specific pre- 
dictions. 

VI. CONCLUSION 

We consider the spin-fermion model for quantum de- 
coherence in solid-state qubits in the pure dephasing 
(i.e. Ti 3> la) situation. We map the evolution of 
the qubit interacting with the fcrmionic environment, 
possibly subject to various 7r-pulse sequences, onto the 
Keldysh path integral. This approach is very general and 
allows one to apply well-developed many-body techniques 
to the problem of the evolution of the qubit coupled to 
the environment and driven by pulses. In the short-time 
limit, we derive the expression for the qubit decoherence 
which involve the product of the noise spectral density 
due to quantum fluctuations of the bath and the filter 
function representing a particular pulse sequence. For a 
non-trivial interacting model of the bath, the Andreev 
fluctuator model, we show that the spectral density has 
1// dependence at low frequencies. Finally, we discuss 
the optimal strategy for the suppression of 1/f charge 
noise by the application of higher-order (beyond spin 
echo) pulse sequences for the problem at hand. One of 
our concrete conclusions of experimental significance is 
that the well-established CPMG pulse sequence should 
be an optimal method for fighting T 2 dephasing when 
the noise spectrum has no sharp ultra-violet cutoff. 
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